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1.  Objective 

Hypergolic  liquids  tend  to  lose  their  self-igniting  capabilities 
at  low  pressure.  One  method  of  counteracting  this  loss  is  to  modify 
the  hardware  in  the  region  where  the  jet  streams  impinge  on  each 
other.  For  example,  when  a  glass  cylinder  one  inch  in  diameter  and 
six  inches  in  length  was  positioned  concentrically  with  the  axis  of 
the  propellant  stream,  the  propellant  was  found  to  ignite  at  pressures 
considerably  lower  than  those  previously  observed  (Ref.  1).  But 
empirical  measures  of  this  type  indicate  neither  the  reasons  for 
improvement  nor  the  extent  to  which  the  performance  can  be  improved. 
Consequently,  they  do  not  revoal  whether  the  desired  performance  is 
attainable  even  under  ideal  conditions-.. 

Instead  of  relying  solely  on  trial  and  error  techniques,  it  is 
really  necessary  to  determine,  from  basic  concepts,  why  hypergolic 
ignition  occurs  less  readily  at  low  pressure  than  at  high  pressure. 


From  such  information  it  may  then  be  possible  to  devise  compensatory 
schemes,  developing  ways  of  assuring  ignition  at  the  lowest  possible 
pressures.  This  report  describes  results  obtained  when  this  type  of 
theoretical  approach  was  applied.  The  actual  procedure  employed 
consisted  of  a  three- step  process.  First,  a  simple  model  of  major 
events  assumed  to  occur  during  ignition  was  developed.  Second, 
equations  describing  these  events  were  established,  and  these  equations 
were  solved  for  the  particular  conditions  treated  by  the  models.  Third, 
the  attempt  was  made  to  identify  the  independent  variables  which  could 
be  controlled  and  which  would  permit  extension  of  the  ignition  range 
to  lower  pressures. 


Before  a  hypergolic  liquid  ignites, a  series  of  events,  each 
lasting  only  milliseconds,  occurs.  This  sequence  is  shown  schematically 
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in  Fig.  1.  Impingement,  atomization,  vaporization,  and  chemical 
reaction  are  shown  as  the  major  events  leading  to  ignition.  Based  on 
this  breakdown,  three  different  models  were  developed,  and  each  was 
individually  subjected  to  further  analysis  in  accordance  with  the  process 
described  above. 

2.  The  Kinetics  Model 

In  earlier  work  (Ref.  1)  the  attempt  was  made  to  show  that  the 
total  energy  released  at  low  pressure  is  less  than  at  high  pressure,  in 
order  to  prove  that  the  process  is  self-sustaining  only  at  high  pressure. 
Yet  equilibrium  thermochemical  calculations  showed  that  the  process 
should  sustain  itself  equally  well  at  all  pressures.  Also,  a  variety 
of  chemical  techniques  were  applied  to  analyze  the  condensed  products 
found  in  the  chamber  after  ignition  was  attempted  (Ref.  2  and  3).  Here, 
too,  chemicals  formed  at  low  pressure  did  not  appear  significantly 
different  from  those  at  high  pressure. 

While  end  products  are  significant  in  affecting  the  ignition 
process,  it  has  become  clear  that  the  rate  of  formation  of  end  products 
is  even  more  significant.  For  the  rate  at  which  the  liquid  decomposes, 
compared  with  the  rate  at  which  it  vaporizes,  determines  whether  the 
decomposition  proceeds  rapidly  or  sk- if  at  all.  Both  decomposition 
and  vaporization  are  rate  processes,  there  is  activation  energy 
associated  with  each,  and  the  two  processes  compete  against  each  other 
for  activation  energy.  Since  the  rate  of  vaporization  of  a  liquid 
depends  upon  the  ambient  pressure,  this  model  implicitly  incorporates 


pressure  as  a  variable.  Based  on  this  model,  the  attempt  is  then  made  to 
calculate  time-temperature  relationships  in  a  droplet  before  ignition 
occurs . 

But  these  time-temperature  calculations  cannot  be  readily  per¬ 
formed.  Values  of  major  parameters  in  the  rate  equations,  particularly 
activation  energies,  are  virtually  unknown.  The  only  data  available  is 
that  which  states  that  some  particular  reactant  liquids  ignited  (or 
failed  to  ignite)  when  mixed  in  the  laboratory  under  certain  particular 
laboratory  conditions.  The  approach  taken,  therefore,  was  to  assume 
values  of  constants  in  the  rate  equations;  then,  to  calculate  the 
resulting  time-temperature  patterns;  finally,  to  pick  out  those  calculated 
time-temperature  results  that  match  the  observed  patterns  of  ignition 
(or  non-ignition)  of  the  measured  periods  of  time,  concluding  that  the 
corresponding  constants  were  appropriate  and  could  then  be  used  in 
further  calculations  as  desired.  The  following  propellants  were  treated 
in  this  way:  UDMH/N204,  MtH/N204,  N2H4/N204,  and  N^/HNOj.  Figure  2 
shows  typical  time-temperature  curves  of  droplets  1000  microns  in 
diameter  when  different  amounts  of  activation  energy  are  available.  When 
the  evaporation  rate  is  assumed  to  be  proportional  to  the  second  power 
of  diameter  (rather  than  to  the  second  power  of  mass)  there  is  a  result¬ 
ant  change  in  the  shape  of  the  curve.  This  becomes  evident  when  Fig.  3 
is  compared  with  Fig.  2.  And,  when  droplets  50  microns  in  diameter 
rather  than  1000  microns  in  diameter  are  considered  there  are  marked 
differences  in  activation  energy  requirements.  This  can  be  seen  by 
comparing  Pig.  4  with  Fig.  2.  These  calculations  suggest  that  large 
particles  do  not  necessarily  show  the  same  ignitibility  as  small 
particles  do.  If  this  ±s  uuu  ignition  ft  hype’-goUc  liquid  at 
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low  pressure  can  be  encouraged  by  forming  the  liquid  particles  in  those 
sizes  which  are  most  ignitible.  Therefore  particle  size  should  be  a 
major  factor  affecting  the  design  of  equipment  upstream  of  the  point  of 
impingement.  These  are  the  main  results  reached  at  this  point  in  the 
investigation,  and  have  been  reported  both  in  Synoptic  form  (Ref.  4)  and 
as  a  complete  manuscript  (Ref.  5). 

The  rate  of  vaporization  of  a  liquid  droplet  is  proportional  to 
some  power  of  particle  size,  but  the  value  of  this  exponent  is  not 
clearly  established.  It  was  necessary  for  these  calculations  (Refs.  4 
and  5)  to  assign  some  value  to  this  exponent.  Based  on  the  particular 
value  assigned  to  this  exponent,  ignitibility  appears  to  vary  inversely 
with  particle  size.  When  a  more  appropriate  value  for  this  exponent  is 
used,  this  conclusion  perhaps  may  be  confirmed.  Yet,  it  is  conceivable 
that  the  reverse  condition  may,  in  fact,  be  applicable,  or  that  the  most 
ignitible  droplets  are  those  in  a  certain  narrow  size  range.  Until 
calculations  with  the  “correct"  exponent  are  performed,  the  relationship 
between  size  of  particle  and  ignitibility  will  continue  to  remain  un¬ 
clear. 

The  appropriate  value  of  this  exponent  can  be  determined  by 
laboratory  measurement.  A  method  of  evaluating  this  exponent  for  any 
liquid  was  developed,  and  a  report  describing  this  method  has  been 
written  (Ref.  6).  The  evaluation  of  rate  constants  of  evaporation  is 
also  treated  in  this  report.  In  Fig.  5,  the  measured  diameter  history 
of  an  evaporating  droplet  of  hydrazine  is  shown.  It  is  clear  from 
Fig.  S  that  a  non-linear  relationship  appears.  But  the  measured  values 
can  be  plotted  in  a  different  way,  with  diameter  raised  to  some  power 
(in  this  case  the  second  power)  so  that  a  straight  line  results,  as 


shown  in  Fig.  6.  When  a  straight  line  is  obtained,  it  is  then  possible 
to  deduce  the  exponent  involved  which  relates  evaporation  rate  with 
particle  size. 

The  rate  of  vaporization  of  a  spherical  droplet  is  shown  by  the 
following  differential  equation: 


d(N/N  ) 

— Z~ 


=  -  k 


(JL f 
j 
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where  N  refers  to  the  mass  of  the  droplet,  k  is  the  rate  constant,  D  is 
the  diameter,  the  subscript  o  refers  to  initial  values,  and  a  is  an 
exponent  which  makes  the  equation  appropriate  for  each  set  of  measured 
data.  Figure  7  shows  curves  of  diameter  versus  time  for  different 
values  of  the  exponent  a.  A  straight  line  is  obtained,  in  a  plot  of 
diameter  versus  time,  only  when  the  exponent  a  is  2.  However,  straight- 
line  relationships  can  be  obtained  when  the  diameter,  raised  to  the 
power  whose  value  is  (3-a)  is  plotted  against  time,  although  when  a  is 
3  then  a  straight-line  relationship  is  obtained  only  if  the  logarithmic 
function  of  diameter  is  plotted  against  time.  These  are  shown  in 
Fig,  8,  9,  and  10,  for  values  of  the  exponent  a  in  the  range  between  0 
and  6. 


In  practice,  both  the  value  of  tho  exponent  a  and  the  rate  constant 
k  can  be  determined  from  measured  values  of  diameter,  and  time,  of  an 
evaporating  droplet.  This  is  done  by  proceeding  in  the  reverse  direction, 
namely,  plotting  diameter,  raised  to  some  power,  against  time,  where  the 
particular  power  that  yields  a  straight  line  relationship  serves  as  the 
indicator  of  the  value  of  the  exponent  a  and  of  the  rate  constant  k. 
Further  details  of  the  method  are  presented  in  Ref.  6. 
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3.  The  Temperature-Distribution  Model 

As  part  of  the  effort  to  improve  the  ignitibility  of  hypergolic 
liquids  at  low  pressure,  this  investigation  turned  to  calculating  tem¬ 
peratures  in  a  propellant  droplet.  In  this  model  there  is  heat 
generated  within  the  propellant  droplet,  there  is  heat  conduction 
within  the  droplet,  and  there  is  loss  of  mass  from  the  droplet  as 
liquid  vaporizes  from  the  surface.  Because  heat  is  assumed  to  be 
generated  uniformly  throughout  such  a  droplet,  while  vapors  are 
released  only  at  the  surface,  a  temperature  gradient  exists  in  the 
droplet  with  the  highest  temperature  at  the  center.  It  is  conceivable 
that  the  temperature  at  the  center  of  a  droplet  can  be  significantly 
higher  than  at  the  surface  so  that  liquid  at  the  center  is  more  ignitible 
than  is  liquid  at  the  surface.  The  difference  in  ignitibility  cannot 
be  estimated,  however,  until  the  difference  in  temperature  can  be 
determined. 

The  heat  flow  equation  for  such  a  droplet,  expressed  in  spherical 
coordinates,  is: 

4^(r2H)  cpH  0  1*  0 

r  r 

where  k  is  the  thermal  conductivity,  ?  is  the  radius,  T  is  the  tem¬ 
perature,  Q  is  the  rate  of  heat  generated  per  unit  volume,  p  is  the 
density,  c  is  the  specific  heat,  t  is  the  time,  and  a  is  the  instant a- 

a 

necus  outside  radius. 

The  droplet  model  assumes  a  constant  surface  t«*3>er*turo,  sc  hat 


the  heat  lost  by  evaporation  is  equal  to  the  heat  conducted  to  the 
surface  from  within  the  droplet. 


da.  .31 
pi  Er 


r  *  a,  t  >  0 


where  X  is  the  latent  heat  of  vaporization  and  where  the  surface  tem¬ 
perature  is  assumed  to  remain  constant. 

A  search  of  the  literature  for  a  solution  to  the  heat  flow  equa¬ 
tion  yielded  only  Awberry' s  solution  (1927).  He  considered  a  droplet 
in  which  heat  is  being  generated,  but  where  there  is  no  evaporation  of 
liquid.  Awberry  reported  the  radial  temperature  distribution  to  be: 
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and  where  T  is  the  surface  temperature,  G  is  the  heat  generation  term, 
a  is  the  outside  radius,  r  is  the  radial  position,  a  is  the  thermal  dif- 
fusivity,  t  is  the  time  and  is  the  initial  temperature.  Awberry  did 
not  even  consider  vaporization,  and  therefore  did  not  take  into  account 
any  change  of  radius  due  to  vaporisation. 

There  foro,  a  method  v;as  developed  based  on  finite-difference 
techniques  to  calculate  not  only  time-temperature  relationships  but  alst 
evaporation-rate  information.  The  model  used  in  these  calculations  is 
shown  in  Fig,  II.  Before  the  equations  were  transposed  into  finite- 
difference  form,  tv*y  were  converted  to  a  non-dimensional ited  fora.  It 
this  way  non-dis&ftuoaal  solutions,  which  are  of  more  general  utility, 
could  be  obtained. 

The  equations  were  non-dimensional ised  by  defining  the  following 
quantities: 
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where  is  the  outside  radius  initially.  After  these  substitutions  arc 
made,  the  heat  transfer  rate  within  the  droplet  is  expressed  by: 


A  further  substitution  is  £ 
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This  leads  to 
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At  the  surface,  the  following  applies: 


da*  r  8T*  ,  ' 
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In  Figs.  12-  13,  and  14  temperature  histories  of  three  points  in 
'he  droplet  are  ; lotted  along  with  the  radius  history  of  the  droplet. 
Figure  12  applies  when  C  =  .01,  Fig.  13  when  C  ®»  .1,  and  Fig.  14  when. 
C  *  1.0.  These  apply  when  the  surface  teaperature  is  the  sane  as  the 
initial  teaperature.  Figure  IS  presents  tice-tetsperature  and  tioe- 


rauius  data  when  the  surface  teaperature  is  different  free  the  ..rdtiai 


teaperature  where  the  difference  is  expressed  by  the  quantity 

ViA)  -  VsA)' 

A  paper  describing  this  work  will  be  subsetted  for  publication. 
Early  work  on  this  subject  wet  s* «d  in  Ref.  ?. 

4,  The  Pre^sure^Dlstribntion 

Also  affecting  the  ignit ability  of  a  propellant  is  the  local 
pressure  is  the  propellant  liquid  after  the  stmas  of  liquid  impinge 
on  each  other.  High  pressures ,  eves  if  local,  should  increase  degree 
of  nixing  of  the  reactant  liquids  and  the  sgoitifcUity  of  the  propellant, 
&  Zi.it  step  toward  shoeing  a  relationship  between  pressure 
distribution  and  propellant  ignitibUity^  a  simplified  aocel  was  assured. 
TWO  circular  jets  opposing  each  other  and  having  the  sane  velocity 
were  considered.  The  flow  was  assured  to  be  iaviscid,  incospressIMe 
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and  irrotational.  First,  streams  of  equal  diameter,  velocity  and  dynamic 
pressure  were  treated.  A  digital  computer  program  was  written,  using 
a  finite-difference  method  to  solve  for  the  flow  streamlines,  jet  surface 
contour,  and  interface  pressure  distribution  (Ref.  8).  The  jet  surface 
contour  and  flow  streamlines  calculated  by  this  program  are  shown  in 
Fig.  16.  Second,  two  jets  of  equal  velocity  and  dynamic  pressure,  but 
unequal  diameter,  were  considered,  as  shown  in  Fig.  17.  Streamlines  at 
various  locations  were  then  calculated  in  Fig.  18.  Interface  pressure 
distribution  and  velocities  are  shown  in  Fig.  IS  for  both  cases  where 
the  jets  are  equal  and,  also,  where  the  jets  are  unequal.  A  detailed 
description  of  the  computer  program  is  available  at  the  University  of 
Santa  Clara. 

After  solutions  for  the  two  axisymmetric  cases  were  obtained,  a 
new  computer  program  was  begun  to  calculate  the  flow  pattern  of  two 
equal-diameter,  equal-velocity  circular  jets  impinging  at  angles  other 
than  180  degrees.  This  is  a  more  difficult  problem  because  it  involves 
a  third  dimension  with  an  additional  degree  of  freedom.  An  initial 
surface  contour  is  assumed,  and  the  velocity  potential  values  on  that 
surface  are  calculated  by  applying  a  known  boundary  condition.  The 
interior  values  of  velocity  potential  are  then  found  by  a  numerical 
relaxation  method,  using  the  Laplace  equation  in  finite-difference  form. 

A  second  surface  boundary  condition  is  then  applied  to  determine  whether 
the  assumed  jet  surface  contour  is  correct.  If  the  boundary  condition 
is  no  satisfied,  the  computer  program  modifies  the  assumed  surface 
contour  until  both  boundary  conditions  are  satisfied.  The  solution  is 
then  complete  for  the  velocity  potential  and  surface  contour.  Gradients 
of  the  velocity  potential  can  be  takan  along  the  impingement  surface  to 
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determine  the  velocity  distribution  there  and,  through  thv  Bernoulli 
equation,  the  pressure  distribution  can  bo  determined.  A  preliminary 
case,  that  of  impingement  at  180  degrees,  wan  solved  with  this  three- 
dimensional  computer  program. 

Since  values  of  jet  surface  contours  must  be  supplied  initially 
to  the  program,  an  experimental  apparatus  was  built  to  measure  surface 
contours.  Unfortunately,  the  program  could  not  be  completed  before 
termination  of  the  contract.  If  this  work  work  were  completed,  it  would 
contribute  significantly  to  solutions  of  three-dimensional  flow  problems. 

There  have  been  requests  for  reprints, for  more  details  about  the 
work  done,  for  a  copy  of  the  axisymmetric  computer  program,  and  for 
information  on  the  status  of  the  three-dimensional  program. 

5.  Conclusions 

Efforts  in  this  project  have  been  devoted  to  improving  the 
ignit -  i:  r;y  of  hypergolic  liquid  propellants  at  low  pressures.  In 
earlier  work,  trial  and  error  methods  were  shown  to  be  remarkably 
effective.  In  subsequent  work,  the  same  goal  was  pursued,  but  the  effort 
consisted  of  work  of  a  more  theoretical  nature,  attempting  to  understand 
the  reason  for  the  loss  of  ignitibility  at  low  pressure  and  for  the 
improvements  obtained  by  empirical  methods.  These  led  to  studies  in 
which  certain  events  in  the  ignition  phenomena  were  simulated.  In  one 
simulation  study,  only  vaporization  and  chemical  decomposition  were 
assumed  to  occur,  and  in  another  study  it  was  assumed  that  only  vapori¬ 
zation,  heat  generation,  and  heat  conduction  occur.  Both  efforts  have 
pointed  to  the  recommendation  that  more  attention  be  paid  to  particle 
size  as  a  factor  affecting  ignitibility.  Another  study  indicated  that 
further  attention  also  be  paid  to  the  angle  of  impact  and  the  pressure 
of  the  impinging  jet  streams. 
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Further  improvements  in  ignitibility  of  hypergolic  liquids  at  low 
pressure  are  destined  to  be  achieved.  It  is  not  clear,  however,  whether 
they  will  be  achieved  through  the  recommendations  presented  here  or  even 
through  the  models  proposed  in  this  research  effort.  But  the  approach 
attempted  here,  of  relating  experimental  observations  with  fundamental 
aspects  of  the  ignition  phenomenon,  should  lead  to  significant  improve¬ 
ments  in  performance  of  hypergolic  liquid  propellants  at  low  pressures. 
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FIGURES 

1.  Sequence  of  Events  in  Hypergolic  Ignition 

2 

2.  Time-Temperature  Relations,  Rate  Mass  ,  1000  microns 
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3.  Time-Temperature  Relationships,  Rate  Diameter  ,  1000  microns 
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4.  Time-Temperature  Relationships,  Rate  Mass  ,  SO  microns 

5.  Measured  Values  of  Diameter  vs.  Time,  Hydrazine 
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6.  Measured  Values  of  Diameter  vs.  Time,  Hydrazine 

7.  Diameter  vs.  Time  for  Various  Alpha  Values 

8.  Straight-line  Plots,  Diameter 3-0  vs.  Time  when  a  *  0,  1,  2 

9.  Straight-line  Plots,  Log  Diameter  vs.  Time  when  a  *  3 
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10.  Straight-line  Plots,  Diameter3”0  vs.  Time  wher  a  *  4,  5,  6 

11.  Finite-Difference  Model  of  Droplet 

12.  Temperature  Histories  (for  r/a,  =0,  0.5,  0.8,  1.0)  and  Radius  History 
when  C  *  .01 

13.  Temperature  Histories  (for  r/a.  *  0,  0.5,  0.8,  1,0)  and  Radius  History 
when  C  ■  ,1 

14.  Temperature  Histories  (for  r/a,  **  0,  0.5,  0.8,  1.0)  and  Radius  History 

when  C  *  1.0  x 

15.  Center  Temperature  History  (r/a,  »  0)  and  Radius  History  when  C  ■  .01  and 
Ti  -  Ts  -  -  i,  0,  1 

16.  Flow  Streamlines  in  Circular  Jet  Impinging  on  a  Perpendicular  Surface 

17.  Impingement  of  Unequal,  Opposite,  Circular  Jets 

18.  Flow  Streamlines  in  Unequal,  Opposite,  Circular  Jets 

19.  Velocity  and  Pressure  Distribution  at  the  Interface  of  Impinging 
Circular  Jets 
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PIG.  1:  SEQUENCE  OF  EVENTS  IN  HYPERGOLIC  IGNITION 
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FIG.  2:  TIME-TEMPERATURE  RELATIONS,  RATE  MASS2,  1000  MICRONS 
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FIG.  3:  TIME-TEMPERATURE  RELATIONSHIPS,  RATE  DIAMETER2,  1000  MICRONS 
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